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Two unusual time-integral conditional regularity results are presented for the three-dimensional 
Navier-Stokes equations. The ideas are based on L'^'^-norms of the vorticity, denoted by Q,rn{t), and 
particularly on Dm = Sim™ , where Om = 2m/(4m — 3) for m > 1. The first result, more appropriate 
for the unforced case, can be stated simply : if there exists an 1 < m < oo for which the integral 
condition is satisfied {Zm ~ Dm+i/Dm) 

Jo \ C4,m / 

then no singularity can occur on [0, t]. The constant C4.m \ 2 for large m. Secondly, for the 
■ forced case, by imposing a critical lower bound on JgDmdr, no singularity can occur in D,n{t) 

for large initial data. Movement across this critical lower bound shows how solutions can behave 
intermittently, in analogy with a relaxation oscillator. Potential singularities that drive Dm dr 
over this critical value can be ruled out whereas other types cannot. 
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I. INTRODUCTION 

The twin related themes of this paper are firstly the regularity problem for solutions of the three-dimensional 
incompressible Navier-Stokes equations and secondly the intermittent behaviour of these solutions. Traditionally the 
first problem, which still remains tantalizingly open, has lain in the domain of the analyst, whereas the phenomenon 
of intermittency has tended to be more of interest to the physics and engineering fluid dynamics communities. It will 
be demonstrated in this paper that these two issues are intimately related & require simultaneous study. 



J"* ■ A. History 

^' 

■ Formally, a weak solution u{x, t) of the three-dimensional Navier-Stokes equations 

idt+u-V)u^iyAu~Vp + f{x), (I.l) 

with divtt = 0, is called regular if the iJi-norm is continuous [l|. What is commonly referred to as 'conditional 
regularity' can be achieved if it is found necessary to impose assumptions on certain system variables such as the 
. ^ . velocity field. The early work of Prodi 0] , Serrin Q and Ladyzhenskaya Q can be summarized thus 0] : every 
• Leray-Hopf solution of the incompressible three-dimensional Navier-Stokes equations with u g L"^ {{0, T)\L'') is 
$_i ] regular on (0, T] provided 2/?- + 3/s = 1, with s G (3, cx)], or if tt e L°° ((0, T) ; L^) with p > 3. The long-standing 
5^ . case s = 3 was finally settled by von Wahl Q and Giga Q who proved regularity in the space C((0, T] ; i^) : see also 
Kozono and Sohr and Escauriaza, Scrcgin and Sverak In summary, the s = 3 case seems tantalizingly close 
to the bounded case s = 2, but not quite close enough. More recent results exist where conditions are iraposcd on 
either the pressure or on one derivative of the velocity field : see the references in Kukavica and Ziane [l2|, [3| j Zhou 
[l3 |. Cao & Titi [H, Ell, Cao Cao, Qin and Titi (for channel flows) Chen and Gala [ll and the revie w by 
Doering [2^. Results on the direction of vorticity can be found in Constantin and Fefferman |2l| and Vasseur j22l |. 
and those on the use of Besov spaces in Cheskidov and Shvydkoy [23|). 

Finally, in recent work, Biswas and Foias [l^l have considered analyticity properties of Navier-Stokes solutions in 
which they have studied the maximal space analyticity radius associated with a regular solution involvin g G evrey- 
class norms. Intermittency properties have also been studied by Grujic [25| and Dascaliuc and Grujic (26l |27{ using 
methods very different from those employed in this paper. 
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B. Motivation and notation 

In an entirely different tliread of intellectual endeavour, the seminal experimental paper of Batchelor and Townsend 
[28j investigated the phenomenon of intermittency in wind tunnel turbulence by observing that the flatness of their 
signals (the ratio of the fourth order moment to the square of the second order moment) took much higher values than 
those expected for random Gaussian behaviour. They concluded that the vorticity is neither spatially nor temporally 
distributed in an even fashion but undergoes local clustering or spottiness, which is consistent with the appearance 
of spikes in the signals interspersed by longer quiescent periods. This is now considered to be a classic characteristic 
of intermittency. These ideas have been developed and extended in many subsequent exp eriments and computations : 
see the papers by Kuo and Corrsin [2^, Sreenivasan [s^l, Meneveau and Sreenivasan j3l| and the books by Frisch [32| 
and Davidson [33| for further references. 

Most of these discussions have been based around Kolmogorov's statistical theory with the widespread use of 
velocity structure functions to study intermittent behaviour. However, structure functions are not easily translatable 



into results in Navicr-Stokcs analysis (see some of the arguments in Kuksin |34| and Dascaliuc and Grujic [27|). The 
main difficulty lies in translating the special conditions needed to prove regularity listed in m Al into sensible physics 
while conversely making sense of the experimental observations in terms of Navier-Stokes variables. The two threads 
can be merged if the spiky nature of the vorticity field is considered in the context of i^'-norms of the vorticity 
u: = curl u with p = 2m 

n^{t)^ (L-^ J \uj\''"'dVj +W0, m>l, (1.2) 

on a periodic box [0, L]^. The additive frequency = vL^^ is present for technical reasons. Clearly the ^m{i) are 
ordered for all t such that 

Wo < fli(i) < fi2(i) < . . . < a„(i) < a„+i(i) . . . , (1.3) 

where L'^/'^ni{t) is the i?i-norm. Control from above on any one of the Vim will also control the ffi-norm from above 
which is the ultimate key to regularity. 

The Navier-Stokes equations have a well-known invariance property under the transformations x' = tx, t' = t^t, 
u — tu' and p = e^p'. Under these transformations fi^ scales as 

, 2777 

f]^" = el7^'" , am = ■ (1.4) 



An(t) - K~'^^™(i)]"'" . (1.5) 



Thus it is natural to define 



While the norms ilm are ordered with increasing m as in (jL3|), the decrease with m. Thus there is no natural 
ordering among the D^- Note also that Di is the square of the i/i-norm. 

For the forced case, the dimensionless Grashof number Gr is based on the bounded-ness of the root-mean-square 
/^„j^ = i^^ll/lli of the divergence-free forcing f{x) and is defined as 

Gr - . (1.6) 



C. Two fundamental results 

There are two results that form the basis of those given in later sections. The first is a theorem on time integrals 
or averages of Dm [13 which uses a result of Foias, Guillopc and Temam [s^. This proof will not be repeated : 

Theorem 1 For 1 < m < oo, and am defined as am — ^m-s ' '"'saA; solutions obey 

f DmdT <c{tGr^ , (1.7) 

Jo 

where rji = Lv^^Eq and Eq is the initial energy. In the unforced case the right hand side is just crji. 
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In [3^ this result was converted into a set of length scales. Let the time average up to time T be defined by 



1 r 

(£'m)T = 1™ sup — / D,n{T)dT (I.. 
D(O) ^ Jo 



D{0) ^ JO 

in which case ()L7[) can be re-expressed as 

{D.m)T<cGr^ +0{T-^) . (1.9) 
Then, motivated by the definition of the Kolmogorov length for m = 1, a set of length scales can be defined thus : 

(LA,-/)'"- := (An)T • (I-IO) 

In [35} the bounds in (jl.9|) were expressed as Re^ instead of Gr^ based on a device of Doering and Foias [stI] who used 
the square of the averaged velocity Uq = (^WuW^^rp to define the Reynolds number Re = UqLi^^^ ; for Navier- 
Stokes solutions this leads to the inequality Gr < cRe^. However, the preference in this paper is to remain with the 
Grashof number Gr. In terms of Re equation (jl.lOp becomes 

L\-^<cRe'^. (1.11) 

When TO = 1, ai = 2, and thus LX^^ < cRe^/^, which is consistent with Kolmorgorov's statistical theory [s^ . 

The second result is a differential inequality for the Dm- Any attempt to time-differentiate the vorticity field 
creates problems because only weak solutions exist. Circumvention of this difficulty requires a contradiction strategy 
commonly used in geometric analysis: assume that there is a maximal interval of existence and uniqueness [0, T*) 
which, for the three-dimensional Navier-Stokes equations, implies that Hi(T*) = 00. In any subsequent calculation, 
if the iJi-norm were to turn out bounded in the limit t ^ T*, then a contradiction would result and so the interval 
[0, T*) could not be maximal. Moreover, it cannot be zero, so T* would have to be infinite. 

Define three frequencies 

where the constants (n = 1, 2, 3) are algebraically increasing with m. The proof of the following theorem requires 
some variations on a previous result (ssj and is relegated to Appendix [A] The dot represents differentiation with 
respect to time : 

Theorem 2 For 1 < to < oo on [0, T*] the Dm{t) formally satisfy the set of inequalities 

Dm < Dl |-W1,,„ (^^^^y + t^2,m| , (1.13) 

where pm = ^m{Am+ 1). In the forced case there is an additive term vas^mGrDm- 



II. A CONDITIONAL REGULARITY RESULT FOR UNFORCED NAVIER-STOKES 

A. Integration of the Dm inequality 

Theorem [2] leads to the conclusion that solutions come under control pointwise in t provided Dm+i{t) > Cp^^Dm{t), 
where Cp^ = [ci.mC2.m]^^''"' ■ The following lemma shows that a time integral version of this controls solutions: 

Lemma 1 For any value o/ 1 < to < oo and e uniform in the range < e < 2, if the integral condition is satisfied 

/ D^^^^dr > Ce,p^ / Df^dr Ce,p^ = [ci,mC2,mY^''"" (II-l) 

Jo Jo 
then Dm(t) obeys D,n{t) < D,„(0) on the interval [0, t]. 

Remark : The case e = 1 turns (jll.ip into /J D,n+i dr > ci^p,^ /J dr both sides of which are bounded above. 
However, the fact that e can take small values suggests a logarithmic result which appears in the following theorem. 
The proof of lemma [T] is included within its proof. 
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Theorem 3 For any value of 1 < m < oo, if the integral condition is satisfied 



In 



l + Z„ 



dT>0, 



(II.2) 



with Ci^m = [2''"" (1 + ci,„iC2,m)] , then D„i{t) < D„i{0) on the interval [0, t]. 
Remark 1 : This result may serve as an alternative to the Beale-Kato-Majda theorem [soj . 

Remark 2 : The exponent p^^^ pulls C4 m down close to 2 for large m which indicates that there needs to be enough 
intervals of time on which Z„i > 1 for (jll.2p to hold. 

Proof: The proof of Lemma [T] is addressed first. Divide p.l3p by Df~'^ and integrate to obtain 



(e) 



Dm 



Df„ dr — ZD. 



(e) 



where Wnjn = (2 — e)wn,m- Noting that pm > 10/3, a Holder inequality then easily shows that 



< 



■'m+l 
D,n , 

•'m+l 



D" 



DLdr 



dr 



DLdr 



pi.SP can then be re-written as 



SlD'm+idr 



(/o 



DLdr 



2,m 



DLdr 



It is clear that no sign change can occur in the denominator of pi.5|) if pi.l[) holds. 

The proof of Theorem III. 21 is now addressed. Divide p.l3|) by D^ and integrate to obtain 



dr 



(II.3) 



(11.4) 



(II.5) 



where we have used (1 + ZmY"^ < 2''™^^ (1 + Z^^^). Re-arranging & using Jensen's inequality 

J exp F{t) dr > exp Q J F(t) dr^ , 
with F = Pm ln(l + Zni), the RHS of pl.6p can be written as 

I f {[l-fZ^]"" -2'''"-i(l + ci,^C2,m)}dr = i / {exp[p,„ln(l + Z™)] -2'''"-i(H-ci,™C2,m)}dT 

> exp 



(11.6) 



(II.7) 



■ / ln(l + Zm) dr 
Jo 



exp [pm In 



(II., 



and thus no zero can develop if pi.2[) holds, as advertised. 

In both cases if a zero cannot appear in the respective denominators then D,n{t) < _D,„(0). It follows that ftm{t) 
is bounded above and thus so is the i7i-norm (r^i). I 
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III. SECOND INTEGRATION OF THE Dm INEQUALITY 
A. A lower bound on Dmdr 
Inclusion of the forcing in (jl.lSP modifies it to 



1 /An + l 



^2,m \ + ^3,mGrD„ 



where pm = ^m{4m + 1). To proceed, divide by (the case e = of Theoreni[T]) to write (jIII.l[) as 



1 d 
2di 



where 



Xm — Wi ■ 



Dm+1 ^ 

"a7 



A lower bound for X„idT can be estimated thus : 



t rt 

Dm+1 dr = 



D 



771 + 1 



< 



Dm"'" dr 



\ 



D„ 



Dt dr 



Dm dr 



and so 



Xmdr > wimt' 



^1 (/o " 



Drn dr 



n73,,„t Gr . 



(IIII.2|) integrates to 



[Dmitf < 



exp|-2 /q Xmdr I 



[DmiO)] - 2w2,m Jo exp {-2 Xm dr'} dr 



Let us recall that p,„ ~ '^m{Am + 1) and let us also define 



2 (m2 - 1) 



then 



(nil) 



(III.2) 



(IIL3) 



(IIL4) 



(III.5) 



(III.6) 



(IIL7) 



Theorem 4 On the interval [0, t] if there exists a value of m lying in the range 1 < m < oo, with initial data 
[D,n{0)] < CmGr^"^ , for which the integral lies on or above the critical value 



Cm{tGr^^^ 4 

where 

V2 > r^iGr^t"^"-^' and where (1 < A,„ < 4) 

A„i = 4 {6„i{l + Prnlm) - Pmlm] 

then Dm{t) decays exponentially on [0, t\. 



m) < I Dmdr 





with 



1 + 4p„t7m 

4(1 + Pmjm) 



< s„, < 1 



(III.8) 



(III.9) 
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Remark : Sm \ 11/20 for large m so enough slack lies between the upper and lower bounds on Dm dr. 
Proof: It is not difficult to prove that 17™ < fi^^j^^ili for m > 1, from which it is easily found that 



Jo Dm dr \ /g Di dr 

Thus (|III.5P can be re- written as 



(/o 






[j*D,dr 





t 

Xmdr > Wl mt ^ ^ , ^ \ p^l^ ^3,mtGr 



> Crnt (nji, „,Gr'^'" - m^^rnGr) (III. 11) 

having used the assumed lower bound in the theorem and the upper bound of J^Didr. Moreover, to have the 
dissipation greater than forcing requires A,„ > 1 so Sm must lie in the range as in (|III.9p because 1 < < 4. For 
large Gr the negative Gr-term in piL6[) is dropped so the integral in the denominator of (jIII.6|) is estimated as 

^ exp (^-2^ Xmdr'^ dr < [2£„wi,„,]"' Gr"^- (l - exp [-2t:^i,™5„,f Gr^"] ) , (III.12) 

and so the denominator of (jIII.6p satisfies 

Denominator > [/^^(O)]"^ - c2.jnC1.jn (2£,„)"^ Gr"^" (l - exp [~2n7i,„£„t Gr^"] ) . (III.13) 
This can never go negative if [£'„i(0)]~^ > ci^mC2,m (2c„i)^^ Gr^^™, which means D„i{0) < CmGr^^"^ . ■ 



B. A mechanism for intermittency 



A major feature of intermittent flows lies in the strong, spiky, excursions of the vorticity away from averages with 
periods of relative inactivity between the spikes. How do these aperiodic cycles appear in solutions and does the 
critical lower bound imposed as an assumption in Theorem |4] lead to this? Using the average notation equation 
pil.6l) shows that if (-Dm)^ lies above critical then D„i{t) collapses exponentially. Experimentally, signals go through 
cycles of growth and collapse so it is not realistic to expect this critical lower bound to hold for all time. . 

Dm(.t) 



1 




Figure 1 : A cartoon of Dm{t) versus t illustrating the four phases of intermittency. 
The vertical arrows depict the region where there is the potential for needle-like singular behaviour. 
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To understand intcrmittency we turn to Figure 1 and draw the horizontal line at Gr as the critical lower bound 
on (-Dm)f Within this allowed range, Dm{t) will decay exponentially fast. Because integrals must take account of 
history, there will be a delay before {Dm)t decreases below the value above which a zero in the denominator of (|III.6P 
can be prevented (at ti): at this point all constraints are removed and the pointwise solution D„i{t) is free to grow 
rapidly again in the interval ti < t < t^. If the value of this integral drops below critical then it is in this interval that 
the occurrence of singular events (depicted by vertical arrows)) must still formally be considered: see the discussion 
in ijlVI Provided a solution still exists at this point, growth in Dm will be such that, after another delay, it will force 
{Dm) I above critical and the system, with a re-set of initial conditions at i2, is free to go through another cycle akin 
to a relaxation oscillator. 



IV. CONCLUSION 



The results of Theorems III. 2 l and HI are new. More appropriate for the unforced case, Theorem lII.2l can be summarized 
thus : if there sufficient regions for which Z„i > 1 on the t-axis, such that 



t /I 



In dr > Zm^ Dm+l/Dm (IV.l) 

V C4_m / 

then Dm{t) < £'^(0). In fact both of the regimes Zm ^ 1 are physically realistic but numerical experiments may 
suggest which of these two regimes are the most manifest and whether a cross-over from one to the other occurs. 

For the forced case. Theorem 2] can be summarized thus: if 

Dm dr > Cm {t Gr^^'- + m) (IV.2) 

then Dm{t) collapses exponentially. The main question is the meaning of the regime below this critical value 

ft 

Dm dT < Cm (t Gr'^"^ + m) . (IV.3) 





It might be natural to suppose, on intuitive grounds, that singular behaviour would be less likely to occur when the 
upper bound is smaller, yet this behaviour cannot be wholly ruled out. What can be ruled out are potentially singular 
spikes that substantially contribute to the time integral of Dm{t) because they would push it over its critical value, 
thereby forcing exponential collapse in Dm{t)- However, there still remains the possibility of needle-like singular 
spikes (depicted by arrows in Figure 1) that contribute little or nothing to the time integral of Dm{t) in pV.Sp . It 
has been shown that the time-axis can potentially be divided into 'good' and 'bad' intervals, the name of this second 
set implying that no control over solutions has yet been found [sl] . To summarize the argument in reference (ssj , it 
is very easy to show that for an arbitrary set of parameters < iirn < 1 , 



a 


Dm+l 




. Dm _ 



Gr-'D^^_,A - rfr>0. (IV.4) 



Thus there are potentially 'bad' intervals of the t-axis on which 

D 



> Cm'Gr-'D';-^, (IV.5) 



D^ 

but on which no upper bounds have been found. Using the fact that Qm+i > ^m, it follows that on these intervals 



Dm+i > CnGri+^-^+i/"™ . (IV.6) 

For large m reduces to D„l+i > CmGr^^^^"^ . Given that /i,„ could be chosen very small these lower bounds could be 
very large indeed, incidentally too large for the Navier-Stokes equations to remain valid. It is possible that these are 
the root cause of the potential singularities discussed above and labelled by vertical arrows in Figure 1. 

It is also possible to interpret this behaviour informally using the so-called /?-model of Frisch, Sulem and Nelkin [4^ 
who modelled a Richardson cascade by taking a vortex of scale Iq = L which cascades into daughter vortices, each of 
scale i„. The vortex domain halves at each step : thus io/in = 2". The self-similarity dimension d is then introduced 
by considering the number of offspring at each step as 2'^, where d is formally allowed to take non-integer values. In 
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d dimensions the Kolmogorov scaling calculations for velocity, turn-over time and other variables have multiplicative 
factors proportional to {£o/in)^^~'^'^^^ ■ see [32|,|4^. Equating the turn-over and viscous times in the standard manner 
one arrives at {£d is their viscous dissipation length) 

LX;^' = eo/id ~ Re^ . (IV.7) 

This gives the Kolmogorov inverse scale of Re^^^ in a three-dimensional domain (43j. To interpret the meaning of 
PV.SP in the light of p.lOp requires the conversion of (jIV.8[) into Reynolds number notation (Gr Re^) as in [33| 



/ DmdT< Cm (t i?e3*'" + 772) . (IV.l 
Jo 



with corresponding length scales 



35 



LA-i < Re^ . (IV.9) 



This makes dm — 2am<5~^ — 1 where the range of 5m is given in pil.9|) . In the large m limit the range of 5m widens 
to 11/20 < (5m < 1 implying that dm lies in the range < dm < 9/11. However, it is conceivable that the sharp 
result for the lower bound on 5m is 1/2 which alters the range to < dm < 1- The upper bound of, or near, unity 
is consistent with the result of CafFarclli, Kohn and Nirenberg [4lj who showed that the singular set of the three- 
dimensional Navier-Stokes equations in four-dimensional space-time has zero 1-dimensional Hausdorff measure. Thus, 
it is possible that singularities that make no contribute to the integral in pV.8[) may conceivably be related to the 
CKN singular set. Whether or not these are physically realisable is open to question. 
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Appendix A: Proof of Theorem [2l 

The following proof has some important variations on that given in [s^. Consider J™ = \uj\'^^'"'dV such that 

— jm= [ |a;p(™"^^u;-{j^Aa; + u;- Vw + curl/} dV^. (A.l) 



2m 

Bounds on the three constituent parts of (jA.l[) are dealt with in turn, culminating in a differential inequality for Jm- 
a) The Laplacian term: Let (j) = ui'^ = u: ■ uj. Then 

< / 0™-iA (lA dV . (A.2) 



Using the fact that A((/)'") = m {(m - l)0'""2|y0|2 ^ ^m-i^^} we obtain 

/ Acjdy < -l(m - 1) / 1 v<?!)|2 dV^ + J- [ A{(j)"')dV 

Jv 2 Jy 2m Jy 

2(m-l) 



|V(w™)rdF, (A.3) 



V 



having used the Divergence Theorem. Thus we have 
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where Am = w"* and ci^m = /{m — 1) with equahty at m = 1. The negativity of the right hand side of (IA.4I) is 
important. Both ||VAm||2 and H^mlh will appear later in the proof. 

h ) The nonlinear term in hA . 1\) : After a Holder inequality, the second term in (jA.l[) becomes 



u>\'^"'\Vu\dV < 



iVtt 



2(m+l~ 



dV 



< c„ 



|2(,n+l)^^ 



|^|2(m+l)^^ 



1/2 Tl/2 



r J ' J 



(A.5) 



where the inequality ||Vm||p < Cp||u;||j, for p G (1, oo) has been used, which is based on a Riesz transform: note the 
exclusion of the case m = oo where a logarithm of norms of derivatives is necessary [39| - see (42| for remarks on 
L°°- estimates. Together with (|A.2p this makes (|A.ip into 

^jn.<-^ I \V{^n?dV + Cm.jlil,JlI^+ I |a;p(™-i)u;.curl/dy. (A.6) 



T/ie forcing term in iA.l\] : Now we use the smallest scale in the forcing I with ^ = L/27r to estimate the last term 
in (|A.6|) 



/ |c^p(™-i)u;.curl/dV^ < ||u;||2™-i||V/| 



(A.7) 



However, by going up to at least > 3 derivatives in a Sobolev inequality and using our restriction of single-scale 
forcing at /c ~ (with I = L/2tt) it can easily be shown that || V/||2m < c ||/||2^~^^ and so 



/ -curl/dy 
Jv 



< cn^J^-^L^wlGr. 



(A., 



d) A differential inequality for Jm '■ Recalling that A„i — allows us to re-write Jm+i as 



Jm+i — ||j4„ 



|2(m+l)/m 
l2(m+l)/m ■ 



A Gagliardo-Nirenberg inequality yields 



1^ 



m\\2{m+l)/m 1^ Cm ||VAm||2 



|3/2(m+l)|| . ||(2)n-l)/2(m+l) 



which means that 



With the definition 



/ p \ 3/2m 

J,n+1 < Cm ( I V(C.'") |2 dv\ J,l?™-1)/ 



(A.9) 



(A.IO) 



(A.ll) 



Pm = -m{m + 1) 

(the factor of | is different from that in [s^l), (jA.lip can be used to form ilm+i 



(A.12) 



(—3 2( 



(„ + l)\l/2(™+l) 



< Cm / |V(w™)|^dF-f TU^ 



l//3„ 



(t~3 J \ 1/2'" 

(-L Jm ) +ruo 



(2m~l)/2{m+l) 



(A.13) 



which converts to 



Cm\L-' / |V(c.™)|^dV^ + n7r > 



(A.14) 
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This motivates us to re-write (|A.6p as 

1/2 



+ c^,n^ln'^-^Gr. (A.15) 

Converting the into fi™ and using the fact that $7™ > ccjo 

° ' C5,m '"^^ I il„ + C6,mn7oGr ^ (A. 16) 



Using a Holder inequahty on the central term on the right hand side, with the definition of Pm given in (|A.12p . (jA.16[) 
finally becomes 

^^n^) " +C3,™Gr ^ . (A.17) 

Given the definition Dm = (wg^^Om)""" it is found that 

o \ ^™ / n \ 



having used the fact that 

(— ^ V™ = 2 (A.19) 

and pm = Pm/oim+i- This converts (jA.171) to p. 131) of the theorem with m = oo excluded. The constants c„.,„ for 

n = 1, 2, 3 grow algebraically with m. ■ 
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